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New crop cultivars will be required for a changing climate characterised by increased summer drought and 
heat stress in Europe. However, the uncertainty in climate predictions poses a challenge to crop scientists 
and breeders who have limited time and resources and must select the most appropriate traits for 
improvement. Modelling is a powerful tool to quantify future threats to crops and hence identify targets for 
improvement. We have used a wheat simulation model combined with local-scale climate scenarios to 
predict impacts of heat stress and drought on winter wheat in Europe. Despite the lower summer 
precipitation projected for 2050s across Europe, relative yield losses from drought is predicted to be smaller 
in the future, because wheat will mature earlier avoiding severe drought. By contrast, the risk of heat stress 
around flowering will increase, potentially resulting in substantial yield losses for heat sensitive cultivars 
commonly grown in northern Europe. 

Most global climate models (GCMs) from the Coupled Model Intercomparison Project phase 3 (CMIP3) 
multi-model dataset predict a decrease of precipitation during summer in Europe (Fig. SI.l) 1,2 . 
Availability of water for plant growth is a key factor determining plant distribution in natural ecosystems 
and is the most important limiting factor in agricultural systems 3 . GCMs also predict a substantial increase in 
temperature and in the frequency and magnitude of extreme weather events in Europe 4 " 6 . It has been shown that 
heat waves in the UK are predicted to increase in frequency (by an order of magnitude), length and severity (peak 
temperature) by the end of the century 7 . Even isolated incidents of extreme high temperature around a sensitive 
stage of crop development, such as flowering, could reduce grain yield considerably, while a continuous period of 
extreme high temperature could result in almost total yield loss 8 . A modelling study for the main wheat growing 
regions of Australia demonstrated that variations in average growing- season temperatures of 2 °C can cause 
reductions in grain production of up to 50%, most of which can be attributed to increased leaf senescence as a 
result of high temperatures 9 . Another study in the USA demonstrated that nonlinear temperature effects on grain 
yields, which show a steep yield decline after temperature exceeds crop-specific thresholds, could lead to severe 
damages to U.S. corn and soybeans yields under climate change 10 . 

The objective of this study was to quantify the impacts of drought and heat stress on winter wheat under climate 
change and to determine which of these factors is likely to increase vulnerability of wheat in Europe. The 
importance of increased temperature relative to changes in precipitation on crop yields was studied before 11 . 
However, these authors investigated the effect of changes in mean climate using a regression model of the crop 
yield sensitivity, which is not appropriate to study the impact of extreme events. In our study, we used a wheat 
simulation model combined with local- scale climate scenarios based on the CMIP3 multi-model ensemble to 
predict impacts of climate change on winter wheat in Europe. We were able to quantify changes in extreme 
weather events and their impacts on wheat, including probability of heat stress around flowering and yield losses 
resulting from drought. By using climate projections from the CMIP3 multi-model ensemble, we were able to 
quantify uncertainty in our predictions and increase robustness of our conclusions, which makes this study 
different from the results reported previously for the UK, where a single climate model, HadRM3, was used 12 . 

The wheat yield could be limited by the grain number and the grain size, which are established to a large extent 
at the period around anthesis (flowering), a stage in development known to be sensitive to high temperature 
stress 8 . If the crop is unstressed, it establishes the grain number and its potential size at sufficiently large values to 
accommodate biomass produced during grain filling. In this case, the yield is effectively source-limited 13 . The 
grain number and the grain size can be substantially reduced if a cultivar, sensitive to heat stress, is exposed to a 
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short period of high temperature around flowering, limiting the 
capacity of grains to store newly produced biomass. In this case, grain 
growth becomes a sink-limited process. In an experiment on the 
combine effects of C0 2 and temperature on the grain yield 
Mitchell et al. (1993) 14 observed that a temperature of 27°C or higher 
applied mid-way through anthesis could result in a high number of 
sterile grains and considerable yield losses. Wheeler et al (1996) 15 
used a temperature gradient tunnel system to demonstrate that a 
temperature of 30°C or higher prior to anthesis significantly reduced 
the grain number and, subsequently, yield of cv. Hereward. In experi- 
ments in Australia 16 , plants were transferred into controlled rooms 
with high temperatures 7 days after the first anthers appeared, show- 
ing that a temperature of 27°C and above could substantially reduce 
the maximum grain size of several Australian wheat cultivars, result- 
ing in substantial yield losses. 

Results 

We used a wheat simulation model, Sirius 17 , combined with climate 
scenarios to predict impacts of climate change on winter wheat 
across Europe. Local-scale climate scenarios were generated by the 
LARS-WG weather generator 18 and were based on the projections of 
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Figure 1 | (A) Anthesis, and (B) maturity. Mean simulated day of anthesis 
and maturity for the baseline (black rectangles) and for the 2055(A1B) 
scenario (box plots) at nine European sites: TR (Tylstrup, Denmark), WS 
(Warsaw, Poland), WG (Wageningen, the Netherlands), RR (Rothamsted, 
UK), MA (Mannheim, Germany), DC (Debrecen, Hungary), CF 
(Clermont-Ferrand, France), MO (Montagnano, Italy) and SL (Seville, 
Spain) (see Table SL1 for details). Box plots are constructed from 15 
predictions of mean day of anthesis (and maturity) for climate scenarios 
(each consists of 300 yrs of daily weather) based on 15 individual GCMs 
from the CMIP ensemble. Box boundaries indicate the 25 and 
75 -percentiles, the line within the box marks the median, whiskers below 
and above the box indicate the 10 and 90-percentiles. 



fifteen global climate models from the CMIP3 multi-model 
ensemble. We estimated the probability of the maximum daily tem- 
perature exceeding a temperature threshold of 27° C and 30° C at two 
stages, at anthesis and at 5 days after anthesis. For heat sensitive 
cultivars, high temperature at anthesis can reduce the grain number 
and heat stress 5 days after anthesis can reduce the grain size. The 
reduction in the grain number and the grain size can subsequently 
reduce the grain yield. To assess the effect of severe drought on grain 
yield, we computed 95-percentile of a drought stress index (DSI95). 
Drought stress index, DSI, is defined as a proportion of the yield lost 
due to water stress: DSI = 1 -Y WL /Y P , where Y WL and Y P are water- 
limited and potential grain yields. DSI95 is a yield loss due to water 
stress, which are expected to occur once every 20 years on average. 

The simulation was run for nine European sites (Table SI.l and 
Fig SI.l). For each site and each GCM from the CMIP3 ensemble, we 
generated using the LARS-WG weather generator 300 years of daily 
weather representing the baseline scenario corresponding to 1960- 
1990, and the future climate scenario corresponding to 2046-2065 
for the A1B emission scenario, denoted as 2055(A1B). For each site 
we selected a typical winter wheat cultivar, which was calibrated 
previously using field experiments (Table SI.l) 19 " 22 . To make a com- 
parison between sites, we used one soil for all sites, medium loamy 
drift with siliceous stones, with relatively low available water capacity 
(AWC) of 1 3 1 mm. Future values of probability of heat stress around 
flowering and DSI95 were calculated for each of fifteen GCMs indi- 
vidually and then presented as box plots to indicate the uncertainty 
in predictions. 

Fig. 1 shows simulated anthesis and maturity dates (day of the year). 
Increase in temperature accelerates wheat development, bringing for- 
ward anthesis on average by 13.2 days and maturity by 17.5 days. Fig. 2 
shows simulated yields for the baseline and 2055(A1B) scenarios. 
Rising C0 2 concentration ([C0 2 ]) will increase the productivity of 
C 3 crops, including wheat. As a result, the simulated yield for the 
2055(A1B) scenario was on average 14.9% higher than for the baseline 
across sites. Although winter wheat cultivars from northern Europe 
are known to be sensitive to heat stress around anthesis 1415 , the lack of 
sufficient experimental data did not allow us to calibrate cultivar para- 
meters for heat sensitivity. Therefore, we assume that all wheat culti- 
vars in simulations were heat tolerant and yields were not affected 
directly by the heat stress around anthesis. 

Accumulation of biomass is proportional to the radiation inter- 
cepted by the plant. Earlier maturity for the 2055(A1B) climate 
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Figure 2 | Grain yield. Mean simulated grain yield for the baseline (black 
rectangles) and for the 2055 (A1B) climate scenario with no change in 
[C0 2 ] of 338 ppm (considering effect of change in climate only, open box 
plots) and increased [C0 2 ] of 541 ppm (gray box plots) corresponding to 
the A1B emission scenario. All wheat cultivars were assumed to be tolerant 
to heat stress around anthesis. 
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Figure 3 | (A) 95-perccentiles of DSI, and (B) soil water deficit at anthesis. 

For the baseline (black rectangles) and for the 2055(A1B) climate scenarios 
(box plots). 



scenario should, in principle, reduce the potential yield, because the 
amount of radiation, which plant can intercept, is lower. However, in 
water-limited environments, the situation is more complex. There is 
a trade-off in maximizing the yield between increasing duration of 
grain filling and avoiding severe water stress. To illustrate the effect of 
changes in climate only on grain yield, we run simulations for the 
2055(A1B) scenario with [C02] of 338 ppm, the same as we used for 
the baseline scenario. On average, the yields for the 2055(A1B) scen- 
ario decreased by about 2.6% (Fig. 2). For the sites where water stress 
was relatively low, e.g. RR or WA, reductions in yields were larger, 
5.6% and 6.7%. However, at some sites, e.g. WS and CF, wheat yield 
increased, because earlier maturity reduced water stress during 
grain filling. 

Fig. 3 A shows box plots of DSI95. For all sites except WA, the 
medians of predicted DSI95 for fifteen GCMs for the 2050(A1B) 
scenario were lower that the values of DSI95 for the baseline scenario. 
This means that despite a decrease in precipitation during summer in 
northern Europe and during the whole growing season in southern 
Europe (Fig. SI.l), relative yield losses from drought were predicted 
to be smaller in the future than at present even for currently grown 
wheat cultivars. Maturity dates are predicted to be 2-3 weeks earlier 
for the 2055(A1B) scenario compared with 1960-1990 (Fig. IB). 
Because soil water deficit (SWD) increases towards the end of crop 
growth, wheat will avoid the most severe drought stress by maturing 
early. It is interesting to note (Fig. 3B) that SWD at anthesis does not 
vary greatly between southern and northern sites and the median of 
SWD at anthesis is predicted to stay at about the same level in the 
future with one exception, CF, where it increased from 50 to 67 mm. 

Medians of monthly mean maximum temperatures for the 
2050(A1B) scenario are predicted to increase between 1.5 and 
3.5°C depending on the month and the site (Fig. SI.l). Fig. 4 shows 
the probability of two events: first, when the maximum daily tem- 
perature within 3 days of anthesis exceeded a temperature threshold 
of 27°C or 30°C (Fig. 4A and 4C), and second, when the maximum 
temperature exceeded thresholds both within 3 days of anthesis and 
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Figure 4 | Probability of maximum temperature exceeding temperature thresholds. Thresholds of 27°C (A,B) or 30°C (C,D) within 3 days of 
anthesis (A,C) or consecutively with 3 days of anthesis and within 3 days of five days after anthesis (B,D) for the baseline (black rectangles) and for 
the 2055(A1B) climate scenarios (box plots). Sirius was used to calculate anthesis dates for each individual year of baseline and 2055(A1B) climate 
scenarios for 15 GCMs from the CMIP3 ensemble. 
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within 3 days of five days after anthesis (Fig. 4B and 4D). 
Probabilities for both events are predicted to increase for all sites, 
except CF, with larger relative increases for northern sites. Exceeding 
a temperature threshold at anthesis will reduce the grain number for 
heat sensitive cultivars, while exceeding a temperature threshold five 
days after anthesis will reduce the potential grain size. Each of these 
events alone can reduce the grain yield for heat sensitive wheat cul- 
tivars. If these events happen concurrently the yield losses could be 
substantial. Although winter wheat cultivars from northern Europe 
are known to be sensitive to heat stress around anthesis 15 , the lack of 
sufficient experimental data did not allow us to calibrate parameters 
for heat sensitivity for individual cultivars and to assess yield losses 
resulting directly from heat stress around anthesis. More experi- 
mental studies to quantify these effects are needed 23 . 

Discussion 

Wheat is the most important crop in temperate zones, including 
Europe, and is the staple food crop for many millions of humans 
and their livestock. New wheat cultivars will, therefore, be required if 
climate change is not to result in losses of yield and food shortages. 
Our results demonstrate that the impacts of changing climate on 
wheat can be counter- intuitive and that the severity of the impact 
will depend on cultivar characteristics and on the spatial and tem- 
poral patterns of climate change. Drought is the most significant 
environmental stress in agriculture worldwide and improving yields 
in water-limited environments is a major goal of plant breeding 24 . 
Some researchers suggest that the impact of drought will increase 
with climate change 25 , emphasising the importance of breeding for 
drought tolerant crops, and this will certainly be true for many crops 
and environments. However, our analysis demonstrates that drought 
will not increase vulnerability of wheat in Europe under climate 
change and relative yield losses from water stress are likely to 
decrease across Europe. In fact, it showed that a more serious threat 
for wheat production in Europe may result from an increase in fre- 
quency and magnitude of heat stress around flowering with potentially 
significant yield losses for heat sensitive wheat cultivars commonly 
grown in northern Europe. Having limited time and resources, crop 
scientists and breeders must select the most appropriate traits for crop 
improvement and should, therefore, focus on the development of 
wheat varieties, which are resistant to high temperature around flower- 
ing, rather than on developing varieties tolerant to drought, which may 
be required for other parts of the world. 

Methods 

To construct local-scale daily climate scenarios, we used the CMIP3 multi-model 
ensemble of climate projections 1 ' 26 . This ensemble was constructed by running several 
GCMs for a common set of climate experiments to assess the uncertainty in climate 
predictions resulting from structural differences in climate models. The direct use of 
climate predictions from the CMIP3 ensemble in conjunction with a crop simulation 
model is not possible because of the coarse spatial resolution of GCMs resulting in 
significant errors and large uncertainty in their output at a local scale, particularly for 
precipitation. The source of errors is related to the fact that many small-scale pro- 
cesses cannot be represented explicitly in GCM, and must be approximated. Various 
downscaling techniques have been developed to underpin studies on regional and 
local scale climate impact assessments, including dynamic downscaling 27 , statistical 
downscaling 28 29 and weather generators 30 . We used the LARS-WG weather generator 
to downscale GCM projections from the CMIP3 ensemble to a local scale 31 . LARS- 
WG has been tested across Europe and showed a good performance in simulating 
diverse European climates 18 . By altering the baseline site parameters using changes in 
climate predicted by GCM, it is possible to generate synthetic daily weather for the 
future. The use of WG in climate change studies allows exploration of the effect of 
changes in mean climate as well as changes in climatic variability and extreme events 7 . 
The latter could be critically important for analysis of agricultural and biological 
systems that incorporate non-linear interactions between system components and the 
surrounding environment 32 ' 33 . Climate projections for each GCM from the CMIP3 
ensemble were treated as equally probable predictions of the evolution of the climate 
system. Only one SRES emission scenario, A IB, and one time period, 2046-2065, 
were used in the simulations. 

The Sirius wheat simulation model was used to calculate crop growth and grain 
yield 17 . Sirius was calibrated for several modern wheat cultivars and was able to 
simulate accurately crop growth in a wide range of environments, including Europe, 



USA, NZ and Australia, and under climate change 20 ' 34 . The model requires daily 
weather data, a soil physical description and management information to calculate 
biomass accumulation from intercepted photo synthetically active radiation day by 
day. Grain growth is calculated using simple partitioning rules. Phenological 
development is linked to the mainstem leaf appearance rate (determined by tem- 
perature) and the final leaf number, determined by responses to daylength and 
vernalisation 35 . 

In Sirius, radiation use efficiency (RUE) is proportional to [C02] with increase of 
30% for doubling in [C02] compared with the baseline [C02] of 338 ppm. A similar 
response was used by other wheat simulation models, e.g. CERES 19 and EPIC 36 . Long 
et al (2006) 37 argued that the results from FACE experiments demonstrated lower 
effect of elevated [C02] on wheat yield 38 ' 39 . To illustrate uncertainty in response of 
RUE to elevated [C0 2 ], we simulated future yields also for a lower value of 15% of 
increase in RUE for doubling in [C02] (Fig. SI3). 

To demonstrate the effect of soil with high AWC on DSI, we run additional 
simulations for a soil with AWC =196 mm (Fig. SI.4). As expected, DSI95 for this 
soil was much lower than for the soil with AWC= 131 mm and showed even larger 
relative decrease for the 2055(A1B) scenario compare with the baseline (Fig. SI.4A). 

Variation in the sowing date could be considered as an adaptation option to a 
changing climate in order to reduce the effect of heat stress or drought on crop yield. 
Fig. SI5 demonstrates the effect of variation in the sowing date on DSI, probability of 
heat stress at anthesis and grain yield at two sites, TR and WS. Sowing dates were 
varied between 262 and 322. The sowing date of 292 for the 2050(A1B) scenario 
produced near-optimal performance of wheat, with the high yield, low risk of heat 
stress and low DSI (Fig. SI5). This date was used as default for the baseline scenario at 
both sites. 
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